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I. INTRODUCTION 



A two-level system (2LS) interacting with a quantum field - electromagnetic field (EMF) in particular - has proven 
to be a very useful model for a wide range of problems from atomic-optical Jl|-0] and condensed matter |||| processes 
to quantum computation JlO[ . For the latter application stringent limits in maintaining the coherence of the the 2LS 
(called qubits) are required . This prompted us to revisit the theoretical structure of the 2LS model, paying special 
attention to its coherence properties. Treatment of spontaneous emission and relaxation are standard textbook topics, 
whereas decoherence and dissipation, especially in the context of quantum computation, are the focus of more recent 
investigations Jll]-[l||]. 

Because of the familiarity of the model (see II. A) and its theoretical and practical values, we do not need to 
emphasize the general motivation, but can go right to the point about the aim and results of this paper. The 
description of this system generally comprises of two parts: 1) Spontaneous emission in the 2LS, and 2) Decoherence 
due to the interaction of 2LS with the EM field, treated as a bath. The first part allows little room for disagreement, 
as it can be obtained from elementary calculations. The second part on decoherence is more subtle. 

Environment-induced decoherenc e |19| has been studied extensively in recent years primarily based on models of 
quantum Brownian motion (QBM) ]2C-|26| for the interaction of a simple harmonic oscillator (Brownian particle) with 
a harmonic oscillator bath (HOB) at a finite temperature, leading to a reasonably good understanding of its charac- 
teristic features. Decoherence of a 2LS in an EM field has been studied by a number of authors, notably |Il] , |i^ , |l6[ , 
and their dissipative and decoherent behavior are reported to be similar to that of a QBM in a harmonic oscillator 
bath. The progression in three stages - quiescent, vacuum fluctuation-dominated and thermal fluctuation-dominated, 
separated by the cutoff frequency and the thermal de Broglie frequency (wavelength) - are indeed characteristic of 
the QBM results |l],E||f5]|2|H . 

Our findings, in contrast, are in stark disagreement from that reported in the literature. We work with the standard 
2LS-EMF model [jl| and obtain an exact master equation for depicting non-Markovian dynamics. Solution of this 
equation for the reduced density matrix of the 2LS shows that the decoherence rate is close to the relaxation rate. 
This is in first appearance rather couter-intuitive, and different from all previous findings. Upon careful deliberation 
we realize that the 'intuition' researchers (including us at the start) have acquired for dissipation and decoherence are 
based on the QBM model which influenced the choice of model in the investigation of decoherence for a 2LS. However, 
we find that such a commonly invoked intuition for QBM in a HOB fails to apply to that of a two-level atom (2LA) 
interacting with an electromagnetic field (EMF) with the commonly assumed type of resonance coupling in quantum 
optics. 



Decoherence in QBM 



Physically, when we say that decoherence of the system of a Brownian oscillator proceeds in a very short time as it 
is brought in contact with an environment, a HOB at some temperture, we are usually conjuring a model with bilinear 
pof (or polynomial [|4|) coupling of the oscillator-bath coordinates, and a ohmic or subohmic spectral function |2l[] in 
the bath. Intuitively, the bath needs to have many degrees of freedom, preferably acting independently of each other 
so incooperatively that the phase information in the system will be dispersed to the largest extent amongst the many 
bath degrees of freedom and affords little chance or takes inordinately long time to be revived or reconstituted (reco- 
herence p9| ) . The opposite picture (of very long decoherence time) is exemplified by two coupled subsystems where 
no coarse-graining is introduced, or for system-environment couplings which maintain some high level of coherence, or 
for an environment whose degrees of freedom have long correlation times like in a zero temperature, supraohmic bath. 
The case of a (spin) particle or (plasma) wave interacting with an averaged (collective) variable from the environment, 
such as the mean field, showing Landau damping in Vlasov dynamics is another example [[30],[5l] • Just as in the spin 
echo phenomena (e.g, Chap. 3 Q]), the basic physics in this case is not dissipation in the Boltzmann sense, but 
statistical mixing p2] ]. We will see that this example is of more physical relevance to our problem than the QBM. 

Coherence in the 2LS 



For the 2LA-EMF system, one clear distinction between an EM field as an environment and a system of harmonic 
oscillators as bath is that the field (coupled to a detector) has an intrinsic spectral density function, not to be chosen 
arbitrarily. For example, it has been shown [^6| that a conformal scalar field in two dimensions coupled to a monopole 
detector has an Ohmic character while in four dimensions it is supraohmic . Barone and Caldeira |33| showed that 
the spectral density function for EM fields with momentum coupling to an oscillator detector is supraohmic. These 
density functions would show very different decoherence behavior from the high temperature Ohmic HOB used in 
many discussions of decoherence, the latter case is what the general folklore is based on. But the most important 
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distinction from QBM is that the 2LA couples with the EMF in the discrete number basis for the field, unlike the 
continuous amplitude basis in the QBM. This fact (which is true in the rotating wave and dipole approximation) 
implies that the 2LS plus EMF system is a resonant one. Hence even though the EM field has just as many (large 
number of) modes as the HOB, only a very small fraction of them in a narrow range of the resonance frequency are 
efficiently coupled to the atom. This is the root cause for the very different qualitative behaviour between the QBM 
and the 2LS as far as decoherence is concerned. 

One extreme case is that of a single mode field described by the Jaynes-Cummings model, where Rabi nutation takes 
place and the atom-field remains largely a coherent system. (For a coherent field, the probability for the atom to be 
found in the excited state at time t regardless of the state of the field) obeys a Poisson distribution. This distribution 
in the photon number induces a spread in the Rabi frequencies, and causes collapse and revival of the Rabi nutation. 
These are distinct features of quantum coherence [Q.) Adding all modes to the field we see spontaneous emission 
and the decay of the atom. The probability of an initially excited atom (remaining in the excited state) decays 
exponentially in the Wigner-Weisskopf form (characteristic of Markovian processes) with relaxation time constant T. 
For purely radiative decay the decay time T\ of the inversion is half the decay time T% of the polarization. There 
is no large order of magnitude differences between dissipation and decoherence time (which in typical QBM high 
temperature conditions could be as high as 40 ]l9|). In fact it is perhaps inappropriate to talk about dissipation for a 
2LA-EMF system because the conditions for a bath to actuate such a process is lacking. The transition from excited 
to ground state is closer in nature to relaxation (in the spin echo sense) than dissipation. In a cavity where excitation 
of the atom from the field (absorption) balances with emission, it is more appropriate to refer to the resonant state 
of the atom-field as a coherent system. In these senarios the distinction between QBM and 2LA cannot be clearer. 



Difference between QBM and 2LS 



So what led earlier authors to make the qualitative claim that 2LS decoheres easily? We think the confusion arises 
when the picture of QBM dissipation and decoherence is grafted on the 2LA-EMF system indiscriminantly. If the 
field which acts as the environment is a phonon field (from ion vibrations, see, e.g., IT^I ) and if the coupling is of 
the non-resonant type, then there is no disagreement. Decoherence should follow the QBM pattern as reported by 
many authors [J Such sources (including atomic collisions in a cavity 0) can be important for some setups. However, 
when one claims that the EM field can decohere a 2LS (with which it is coupled in a resonant way, as in the standard 
model) that is where we disagree. 

Quantitatively, the model for the 2LS used by most authors for the discussion of decoherence inspired by QBM 
type of behavior has the atom in a a z state (the diagonal Pauli matrix) coupled to the field mode operators b* , b. This 
type of coupling term (call it a z type for convenience) commutes with the Hamiltonian of the system, and admits a 
diagonalization in the eigenbasis of the Hamiltonian. The field is coupled to the atom as a whole and thus is insensitive 
to the 2 level transition activity. In particular it does not probe the resonance or coherent properties of the two level 
atom, which is the most important feature, for quantum computation. By contrast the standard model for 2LA-EMF 
which we studied has a a± coupling (call it standard coupling) to the field modes which highlights the 2 level activity 
of the atom and the field. This coupling considered in the standard model is indispensible, i.e., it cannot be removed 
from the two-level atom as it defines it and will be present in any realistic situation. What then is the origin of the 
QBM type of contribution to the 2LS ? 

If one accepts an environment other than the EM field, the question comes down to the characteristics of the 
experimental apparatus. For well- prepared ion traps we would expect it to be rather unimportant. If the EM field 
is the only field present, we can still ask if a QBM type of coupling term with the EM field would appear, and if yes, 
how strong would its effect be? This would be a useful way to accomodate the two different types of coupling terms. 

Recall that the standard model is derived under the dipole and rotating wave approximation. In the next section we 
will show that the o~ z type of coupling appears only in the next order expansion after the dipole approximation. Since 
these are good approximations for a large class of atomic states when the atom is nonrelativistic, the contribution 
from the QBM type of coupling used in [|ll 12 Iff] should be negligible and its ensuing decoherent effect insignificant. 
In this sense the EM field does not in leading order of approximation act like a bath in the QBM way, and coherence 



1 Even in such cases, one also needs to pay closer attention to the QBM behavior than what has been accorded for this model. 
Subtle points unnoticed before include, e.g., the imposition of a high frequency cutoff and Ohmic spectral function which 
restricts to a Markovian behavior J2l| can lead to a violation of the positivity of the reduced density matrix , the violation 
of the fluctuation-dissipation relation J2j| , and the prolongation of coherence in a low temperature supraohmic bath |m| . They 
deserve more attention in the theoretical design of cavity qubit computers. 
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in a 2LA-EMF system is quite well preserved (excepting other processes, e.g. ^SlQ ). 

Our puzzle over the result on decoherence in the 2LS reported in the literature was what prompted us to begin this 
study. Without letting any familiar and convenient analogy influence our judgement, and without any preconceived 
notion, we choose to perform a first-principles calculation of the two-level atom (2LA)-clcctromagnctic field (EMF) 
system making as few assumptions and covering as wide a range of conditions as possible. We use the influence 
functional method pp|| to take into account the full backreaction of the field on the atom, while adopting Grassmannian 
variables for the 2LA and the coherent state representation for the EMF. We find exact master equations for the full 
(non-Markovian) dynamics in the cases of a free quantum field and a cavity field at zero temperature. 

In the next section we present the model and the formalism. A detailed derivation of our model is contained in 
Appendix A. In Sec. 3 we derive the master equations. In Sec. 4 we study different mode composition of the field, 
including that of an atom in a cavity. We end in Sec. 5 with a discussion of our findings and their implications. 

This is the first in a series of papers on 2LAtom and quantum decoherence. The subsequent papers will treat 2LA- 
EMF interaction at finite temperature, for EM fields in a coherent and squeezed state, and for multipolar models 
(where coupling other than the minimal is assumed). We will also tend to collective qubit systems and moving 
atoms interacting with an EM field. These results will have corresponding applications in atom optics and quantum 
computation problems. 



II. THE INFLUENCE FUNCTIONAL 



A. The Model 



Our model for atom-field interaction is the standard one (see Appendix A for details) fl],|],[§ The total Hamiltonian 
for a (stationary) atom interacting with a quantum electromagnetic field (EMF) under the dipole, rotating wave (RW) 
and two-level (2L) approximation is given by 



H = hw S z + hJ2 [^blk + (<?kS+S k + ffkS-S k )] (ii.i) 

k 



where 6 k , 6k are the creation and annihilation operators for the kth normal mode with frequency Wk of the electro- 
magnetic field (thus for the field vacuum 6k|0) = 0, [6k, 6 k ,] = <5k,k', for all k.), and uiq = u>2i is the frequency between 
the two levels. Here 

S z = -zfrzi &± = a± = -(<7 X ± i<T y ) 

where <J x ,y,z are the standard 2x2 Pauli matrices with a z = diag(l, —1), etc. The coupling constant <?k = d2ik/k(X) 
where 

d -^~V&V dir ^ (IL2) 

and dij = e J (j>iX(f)jd x is the dipole matrix element between the eigenfunctions <pi of the electron-field system, ek CT 
is the unit polarization vector ( a = 1,2 are the two polarizations), and /k(x) is the spatial mode functions of the 
vector potential of the electromagnetic field (in free space, /k(x) = e~ akx , V is the volume of space.). Under the 
dipole approximation /k is evaluated at the position of the atom X. Since dy = dji, d^k = <^jik, we will choose a 
mode function representation such that is real. 

To see how this could possibly be related to the a z type of coupling with Hamiltonian (used by e.g., |12 1(| for the 
study of decoherence in 2LS) 



H = huj Q S z + [w k S k S k + ha z (#k6k + 5k6 k ) (II.3) 

k 



2 Our Hamiltonian is given in the so-called minimal coupling (MC) as different from the multipolar coupling (MP) which 
may be more relevant to atoms in a cavity because the expicit Coulomb interaction between the atom and its image charge is 
removed. 
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we examine the next term after the dipole approximation in (A. 15). This has a contribution to g^ even when i = j 
which is equal to 



3iik = c k k • qi 



where 



and Ck is a constant given by 



<l t =^2 J <f>%5x(p ■ e klJ )<j}idx 3 



Ck = (2ftw k e V)" 

m 

This generates an additional coupling term 

^2 VzigikK + Sik^k) + l(.92k^k + .92k^ k ) 

k 

where 

3lk = Silk — 522k, 32k = .9llk + 322k 

This gives the lowest order a z type of coupling in a 2LA -EMF system. The ratio of the coupling gik of the a z type 
in Eq. (II.3)to the dipole coupling in Eq. (II. 1) is 

I , | , k(qi -q 2 ) ^ k |qi -q 2 | m . 

Wffk = I ~, 1 < -, (II.4) 

mCJkdi2 TOWk«12 

Thus the a z type of coupling generated from the 2LA- EMF interaction will be significant only for very high frequen- 
cies u>k of the EM held, a point intuitively clear from the meaning of the dipole approximation. 



B. Grassmannian Variables and Coherent State Integrals 

Since Feynman and Vernon |20[ invented the influence functional method this formalism has been applied to treat 
the Brownian motion of a harmonic oscillator interacting with a harmonic oscillator bath by many authors |2l],|^,|22|,Q . 
The two - level system in tunneling has been discussed in detail by Leggett et al ||, but the derivation of a master 
equation by this method which can traverse the non Markovian regimes has not yet been carried out. We shall perform 
such a calculation for a two level system, with the aid of Grassmaninn variables convenient for treating fermions, and 
the coherent state representation in a path integral form. We construct the coherent state of the combined atom-field 
system as 

|{z},ry) = \{z}) x \r,) (II.5) 

where \z), z a complex number, denotes the EM field coherent states and \rj, r\ a Grassmannian or anticommuting 
number, denotes the electron coherent state. The transition amplitude between the initial state (i) at t = and the 
final state (f) at t = tf is expressed formally as [|4j (here we suppress the index k, 

(fj f , z f ; t\rn, Zil 0) = J DzDzD V Dfje* s ^ z ^ {11.6) 



where the action is 



z It 

-S[z,z,rj,fj} = zz(t) + fjr}(t) - ds(zz + fjf) + ~ H(fj,r),z,z) (II.7) 



Here H is the Q-symbol of the Hamiltonian pBl and there is an implied summation over field modes. 
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H{rj, r], z,z) = h ^2 ^kZkZk - g k (zkV + V z k) + ^oVvj (H.8) 

In (II. 8) we have substracted a constant term ^wol to (II. 1) so that the ground state now has zero energy. Henceforth 
we set h = 1. 

The Hamiltonian in equation (2.8) is not a c-number function; it has terms that are odd. One might then question 
the validity of equation (2.6) for the path integral; it clearly exists as a formal expression, but its evaluation with a 
saddle point method, that is based on the Hamiltonian of equation (2.8) might be problematic. We dispell this doubt 
for the vacuum case with an operaytor method proof of the master equation in Appendix B. It shows that at least for 
the vacuum case tha saddle point evaluation yields the correct result. The general cases need separate considerations. 
For many qubits coupled to the EM field vacuum, we believe that the path integral method yields a simpler treatment 
than the operator method. 

The integration is over all paths satisfying 

z(0) = z t z{t) = z f (II.9) 

v(0)=Vi m=v f (H.10) 

We assume initially that the density matrix of the total system+environment is factorizable p(0) = p e {0) <8> pb{0). 
Only at that time would z and rj be pure complex and Grassmannian numbers respectively. As the system evolves, 
both r] and z contain Grassmann and c-number parts. The mixing of even and odd parts (note gk is odd) comes 
about as the initially factorized atom state becomes " dressed" . 

In the open system philosophy, as we are interested in the averaged effect of the field on the atom, the atom is 
considered as the 'system' while the field as the 'environment'. The path integral is performed over the variables z, 
while r], fj are treated as external sources. When only one field mode is considered, we have 

(zf, t\zf, 0} Vtf j = J DzDzcxp^zz(t) — J ds(zz + i(cuzz — g(zr\ + fjz)(s)^ (H-H) 

with summation over paths satisfying the boundary condition (II. 9 ) for z. We use the saddle point method. Mini- 
mizing the action yields the following equations 

z + iuoz^—igrj (H-12) 
z — iuoz = igfj (11.13) 

with solutions 

z( a ) = z ie - luJS - ig f ds'e-^ s - s '^(s') (11.14) 
Jo 

z( s ) = z f e- luj(t - s) + igj^ ds'e- l ^ s - s ' l fj(s') (11.15) 
Using these for the transition amplitude (11.11 ) with the minimum value for the action, we obtain 
(zf,t\z i ;0) rit fj=exp^ZfZie~ iwt -ig z f J dse' luj{t ^ s) T](s) 

+ [ dse~ lujs f](s)z l -g 2 [ ds [ ds'e-^ s - s '\fj(s)r]( S ') \ (11.16) 
J Jo Jo J 

A prefactor in the coherent state path integral is equal to one. Now the influence functional due to this single mode 
reads 

r[r,,fj;rff?] = [ d ~ ZldZl d ~ Z ' dZ ' % d *l dz l c -™--*M--*"> 

' J IT TT 7T 

x(zf,t\zi;0) Viij (zi\p \z'i) (^jOl^;*)^ (H.17) 
where the completeness relation for (unnormalized) coherent states has been used 

dzidzi 



J 



z \z)(z\ = 1 (11.18) 

TT 
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Writing with an obvious identification 

{zf,t\zi;0) v>f j = exp(AzfZi +izf/3 + ijZi + D) (11.19) 

we can use the identity 

J dz^Zi e _ Ez+f - z+Ef = eff (n 2Q) 

to obtain 

mv,ff,f/] = ^ D+D,) (11.21) 

for an initial vacuum state p = |0)(0|. Substituting, we get the contribution to the influence functional from one 
mode 

r k [fj,V,fj',v'] =expj<7^ ds J d S '[f)'( S ) V ( S ')e-^ (s -^ 

+f]'(s')r 1 (s)e lu ' k{s - s " ) - fj(s)r 1 (s')e- Wkis - s ' ) - fj' (s')r](s)e Wk< - s - s '^ } (11.22) 
The influence functional for all modes T — J\ k Tk is finally given by 

F[fj,r),fj',r)']=exp^j ds J ds' (fj,(s - s')[fj' (s) - fj(s)]r](s') 

+f(s-s')f?(s')[ri(s)-T/{8)])} (11.23) 



in terms of the kernel 



m«) = ( IL24 ) 



k 



III. THE MASTER EQUATION 
A. The reduced density matrix propagator 

Having computed the influence functional we have an expression for the reduced density matrix propagator 
JiVfVf: tlv'iVi'i 0) = J DfjDrjDr]' Drj' exp |^r?(<:) + ffvi'it) — J ds(fjfj + ffrf + iufjrf — icvfj 1 r\ (s) 

+ fo ds [ ds ' {il[s ~ s ' )W[s) + ^ sMs ' ] + M * (s _ 8 'W S 'M*) + »/(*)])} (ni.i) 

where the summation over all paths obey the boundary conditions (11.10) and 

?j'(o) = ff. v '(t) = v ' f (IH.2) 

We can compute the path integral with saddle point evaluation and get 

fj + iur) + [ ds'fj,(s - s')ri(s') = (III.3) 
Jo 

fj'-iufj'+f ds'Li*(s-s')ri'(s') = (III.4) 
Jo 

fi' + iwr]'+[ ds' fi(s - s')ri(s') - f ds' /j,*(s - s')[r)(s') + r]'(s')} = (111.5) 

Jo J s 

fj-iufj + f ds'/j,*(s-s')fj'{s')- [ ds'fj,{s- s')[fj'{s')+fj{s')} = (III.6) 

Jo Js 



It will turn out that only the solution of the first two of these equations will contribute to the path integral. We will 
therefore write 



r)(s) = rjiu(s) r\ (s) = rj' f u(s) (HI.7) 

where u,u are the solutions to equations (III. 3), (III. 4) under the condition 

u(0) = u{t) = 1 (III.8) 

Now equation (III. 3) is a linear integrodifferential equation of first order and as such can be solved with the use of 
the Laplace transform and the convolution theorem. It is easy to show that 

u(s) = £-i ( —— ) = -— / - . (III.9) 

\z + tco + n{z) J 2m J c _ ioo z + IUU + fi(Z) 

where jl(z) is the Laplace transform of the kernel (11.24) and c is a real constant larger than the real part of the 
poles of the integrand. It turns out that this function u(s) contains all necessary information for the computation of 
the density matrix propagator. Substituting our expressions (III. 8) and (III. 9)) in (III. 1 ) we can obtain the following 
expression for the propagator 

JiVfVf-, *fe 0) = exp (f} f T]iu(t) + fj' lV ' f u(t) - [1 - u(t)u(t)]fj' im ) (111.10) 

Since we are using coherent state path integrals we have departed in our evaluation from the standard saddle point 
approximation used on configuration space path integrals. In these cases , the standard procedure is to distinguish 
the imaginary part of the kernel as corresponding to dissipation and consider only its contribution when performing 
the saddle point evaluation. The resulting equations are then the classical dissipative equations of motion. But in the 
case of the coherent state path integral, there is no correspondence between extrema of the action and actual classical 
paths. Hence, there is no sense in splitting the kernel rj into real and imaginary part, and the saddle point evaluation 
should be carried out for the whole of the exponential. 



B. Master equation for a field in a vacuum state 

It is a standard procedure now to find the master equation |2l],|24|]. We compute the time derivative of the propagator: 

j = {ufjfrh + iLfjfrf + -±-±ff irH )j (in.ii) 

The next step is to remove from the above equation the dependence on the initial values. This is done with the use 
of the following identities 

U 01] f U 07] j 

Note that we are suppressing (for ease of notation) symbols denoting left or right Grassmann differentiation. In all 
our expressions we implicitly assume that differentiation with respect to rj is always right and with respect to f) always 
left. 

For the density matrix at time t 



Pt(Vf,Vf) = J dj]idr]ie nm dvi'^e J^/r^; ^77-77,; 0)po{r]i>Vi) (111.13) 
we obtain the evolution equation 

d u 80 uSp 4z(uu) 5 2 p 

at u or/ u 01] uu 07/077 

This is one of our main results: The master equation for the two-level atom interacting with a an environment of 
electromagnetic field at its vacuum state. The effect of the field is contained within the function u which can be 
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determined by the solution of equation (III. 3) or cquivalently by the computation of the contour integral (III. 9). In 
the next section we are going to find explicit expressions for u for particular choices of the field configuration. 

Let us return for the moment to equation (III. 14) and write this in an operator language. It is easy to verify that: 

4H + " m= s+pS - (IIU5) 

where S + = (1 + er z )/2. If we write 

= r(t) + ifi(t) (111.16) 



«(*) 



the master equation reads 



^p = -i[Sl(t)S+S-,p] +T(t){S + S-,p} - 2T(t)S-pS+ (111.17) 

where 

H(t) = fi(t)E+ (111.18) 

The first term corresponds to the unitary Hamiltonian evolution, only now the effect of the environment has induced 
a time dependent shift in the value of the frequency, the second term is time dependent dissipation and the third 
corresponds to noise. 



C. Spontaneous emission 



To show how the standard results are regained, and to understand the meaning of the new function in the master 
equation, let us consider the physical process of spontaneous emission. Start with a generic initial density matrix 

P={ 1 -/ l) (IH.19) 

its corresponding Q-symbol is 

p(Vi V) = x + V* 1 ! + UV + (1 — x)fjr) (III. 20) 

If we evolve it with the density matrix propagator (III. 10) we obtain for the state at time t 

Pt(fj, rj) = 1 - uu(l - x) + {uy*rj) + {ufjy) + (uu(l - x)) fjrj (III. 21) 

corresponding to 

/ uu(l -x) uy \ 
pt ~\ uy* 1 - uu(l -x)j (IIL22) 

Considering the case x = y = we get for the probability of spontaneous emission 

P(l->0,t) = l-uu (111.23) 

Also we should remark that the rate of decoherence in the energy eigenstates is governed by the absolute value of the 
function u (the off- diagonal terms). But on the other hand u itself determines the rate of nergy flow from the atom 
to the environment. Hence for our particular choice of initial state (vacuum) we find that decoherence and relaxation 
time are essentially identical . We shall use this equation to study decoherence in an ion trap in a later paper. 



IV. FIELD MODES AND ANALYTIC U(T) 

Our master equation (III. 17) depends solely on the function u{t) , which in its turn is determined by the kernel 
p(s) . In this section we will try to give some analytic expressions for this function in various different cases. 
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A. A single mode 



To connect with known results Q, let us start with the case when the field contains only a single mode with 
frequency u>k — k. Then /i(s) will read 

H(s)=g 2 e- iks (IV.l) 

and 

r°° n 2 
m = -9j Q ^e^ = JL- k (IV.2) 

The integrand has two poles at the solutions of the equation 

z 2 + i{uj + k)z-uk + g 2 = (IV.3) 

given by 

.LU + k±[(LJ~k) 2 +g 2 }^ 2 . 

z = -i — — — — = -twi,2 (IV- 4 ) 



Hence 



u(s) = k ^ e -^s _ k - u 2 e _ iU2S 

UJ2 —W\ LU 2 — LUl 



This result is in agreement with standard ones 



B. Infinite number of modes 



Now we consider the case of the vacuum electromagnetic field in free space i.e. not constrained by a cavity. The 
kernel will read then (using equation (11.24)) 

/rl 3 k- \ 2 r°° rl 

JL* fc-i e -*« = L. ^ kdke- ik ° = ± V (a) (IV.6) 

where 

i a 2 r°° 

u(s) = / dke- lks (IV.7) 

Note the factor of 2 in (IV.6) coming from the two photon polarisations and that we in view of (II. 3) we have written 
.9k = Aw k 

Since the integral (IV.7) is not convergent, we will introduce an exponential cut-off in the higher frequency modes. 
The presence of the cut-off is of physical significance since we do not expect high electromagnetic modes to couple 
with our two-level atom. 

Hence the kernel v will read 

v{a) = At / dke~ lks - ke = (IV.8) 



n 2 



The Laplace transform of v is then 



\2 foo -sz \2 

v{z) = -j / ds- — - = -— e- lcz Ei{-iez) (IV.9) 



where Ei denotes the exponential integral function analytically continued to the complex domain. At the limit e — ► 
this is essentially 

Ei(-iez) = 7 + log(-iez) + 0(e) (IV.10) 
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where 7 is the Euler-Macheronni constant and the logarithm is taking values in the primary branch. Thus jl(z) reads 
(i^(0) is here v(s = 0) obtained by the integration by parts of the Laplace transform) 



j\2 \2 

fi(z) = -i/(0) + zv{z) = 5 -ze- tez Ei{-iez) (IV.ll) 

Note that the cut-off e affects significantly fJ-(z) only at large values of z, which essentially correspond to the very 
short time limit, i.e the time where the two-level atom starts "getting acquainted" with the photon reservoir. At 
larger times (t » e) we do not expect the cut-off to contribute significantly in the evolution. This is a rather typical 
behaviour in quantum Brownian motion models, provided that the ultraviolet cut - off of the environment is much 
larger than the natural frequencies of the system. 

To evaluate the integral we first have to find the poles of the denominator. We can do that in a perturbation 
expansion . First let us absorb the divergent z/(0) factor in a frequency renormalisation. 

\2 

u> = u - -=- (IV.12) 

so that we need find the zeros of the function z + ituz — ^e~ tez Ei(—iez). Looking for a solution in the vicinity of 
z = —id) we find 

z = —iiu — fh(—iu>) 

-i(u> - ^4- log(e^)) - — + 0(A 4 ) := -in - F (IV.13) 

We have a pole with a negative real part and we can verify numerically (also physically expected) that there is no 
pole with greater real part. This means when evaluating the inverse Laplace transform we can ignore the contribution 
of the branch - cut at z = (being on the right of the pole) and hence after some time, where all possible other poles 
with absolutely larger value of their real part will have stopped contributing the solution will be 

u(s) = e~ ins - rs (IV. 14) 

This implies a Markovian time evolution and the identification of decoherence - relaxation time with = This 
is a general feature of the presence of a continuum of modes as can be seen from the case of an atom transparent to 
all modes but the ones in a strip , say [ui,w 2 ] containing the resonance frequency. 
It is easy to verify that in this case we have again 

^) = -^( W2 -wi) + ^log(^^) (IV.15) 



Hence defining again 



0J\ — iz 



A 2 

u = u ^(cj 2 — <^i) (IV. 16) 

7T 



* = - ^ log f ^^r) )- — + 0(A 4 ) (IV.17) 



we can find the pole at 

•(- X2uJ ^ ( 

Z = —tlld — loe- I 

LU — LU — I J 7T 

Note that the real part of the pole comes from the presence of a minus sign in a logarithm of some real valued object. 
Hence in the case where the atom 's frequency is outside the strip of interacting modes there will be no dissipation. 
This feature separates us from the QBM case, characterising the atom - field system as primarily a resonant one. 



C. Atom in a cavity 



Let us now consider the case of the atom lying within a cavity consisting of two parallel plates at distance L. The 
field satisfies Dirichlet boundary conditions on the surface of the plates, hence the modes in the normal direction to 
the plates are multiples of w/L. The kernel then reads 
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/i(s) 



2nL 



E 



kdk 



o (fc 2 + (nir/L) 2 ) 1 / 2 



-j(fe 2 + (r«i-/L) 2 ) 1 '' 2 (s-ie) 



— E 

2-kL ^ 



dke- lk( - s ~ le) 



t/L\ 



!tt 7", p -I- is -^-^ 



—i\mz/L\ (s— ie) 



27ri e + is 

n 

1 11 + £-"r/L(s-ie) 

~ 2ttL e + is 1 - e -"r/L(i.-ie) 

Hence we can compute its Laplace transform 

, . -iX 2 r°° dse- sz 1 + e -™/L(.s-i') 
fi{z) 



(IV.18) 
(IV.19) 



2-kL j 
-i\ 2 



2nL JHeZ '2L; ) 



it \ g — iir / L(s — it) 

in 



where J(x, a) is denned by 



J(x, a) 



dye~ v 



coth(ay) 



(IV.20) 



(IV.21) 



and appears in equation (IV.20) through analytical continuation in the complex plane. This integral can actually be 
computed at the limit of vanishing a; (e — > 0) - see reference pffl, equation 3.427.4. 



-i\ 2 
ir 2 e 

irL 



X 2 

— -z log(«e 7 e2) 

^,Lz. Lz .Lz. Lz 1 Lz 

io g r — - — io g (— + — + - log — - 2 

ITT lit ITT lit A 2«7T Z 



0(e) 



(IV.22) 



Note that p,(z) is a sum of the term of case 2 (L — > cxd) and a finite one (no dependence on e ) . The logarithm of 
the T-function gives a countable number of branch cuts at z = n positive integer, the resonance modes of the 

cavity . Again the important pole has a negative real part. We can again compute the pole perturbatively. It lies at 



A 2 A 2 
-i I Q H — - H) •\og(e 1 'euo) logr(— Llv/tt) 



Llu _ Lu 1 Lib . 
log {-Lu tx) log (- — 

7T IT 2 Zir z 



(IV.23) 



Clearly the logarithm of the gamma functions is the term out of which the real part of the pole appears. Since the 
real part of the pole is negative the branch cut is excluded from the integration contour and hence 



u(s) 



—iQs—Ts 



(IV.24) 



Here T gives a dissipation constant. . In Figure 1 we give a plot of the real part of the pole (— T) as a function of the 
frequency u>. Note that it has sharp maxima on the resonance points, implying persistence of coherence . 

Already from the approximation (IV.23) we observe that the difference in the renormalized frequency from the 
case L — > oo (Aw) is finite. Unfortunately perturbation expansion is not reliable when Q is close to the resonance 
frequencies (this corresponds to negative integers arguments in the T-function where it diverges) and for this regime 
we have not been able to get any analytic results. In Figure 2 we have plotted the dependence of Aw = Sl[L] — f2[oo] 
as the frequency changes. This effect of the frequency shift for an atom within the cavity is well known, as well as its 
relation to the Casimir effect Ipql. 
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FIG. 1. Real part of pole (negative dissipation constant V) as a function of ui/(n/L). 
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FIG. 2. Frequency shift AS! as a function of u>/(n/L). 



V. DISCUSSION 



Let us now integrate what we have found and look at the overall picture. The physics of a 2LA-EMF system at 
zero temperature is characterized by a number of time constants: 

1) The inverse natural frequency uJq 1 

2) The inverse coupling constant g^ 1 — y/uJk/\ 

3) The relaxation time constant 

4) The cavity size L (divided by c) 

First consider a zero temperature field in free space, thus ignoring factors 4). Start with only one mode in the field 
in resonance with the atom, then the system undergoes Rabi nutation with frequency fl « gyfn 
photon number in the field. The collapse time (assuming a large mean photon number h) is g~ 



+ 1, where n is the 
, and revival time is 
but shorter than g . 



(This is the 



2ny/E/g. Atom excitation becomes significant in a time much greater than i 
condition for a first order perturbation theory to give reasonable results.) For a large number of modes, spontaneous 
emission occurs at the relaxation time scale = n/g >> lUq 1 which we found to be the same as the decoherence 
time - the time for the off diagonal elements of the reduced density matrix to decay (Sec. 4). When the mean 
number of photons in the field is large (n >> 1), they become comparable to the collapse time. This is a measure 
of the coherence in the atom-field system, and is controlled mainly by their coupling and the photon number in the 
field. We see that with the resonance condition, the nature of decoherence in 2LS is very different from the QBM 
situation, where phase information in the Brownian particle is efficiently dispersed in the many modes in the bath 
coupled almost equally to the system. As we remarked in the Introduction, the identification of the phase information 
and energy flow from the 2LS to its environment is similar to the spin echo phenomena (Landau 'damping') which is 
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based on statistical mixing rather than dissipation. The mathematical distinction lies between considering the system 
coupled to the discrete number basis (our model) and the continuous amplitude basis (QBM) of the environment. 
The latter case essentially produces noise that drives the system in a way insensitive to its own intrinsic dynamics. 
While in the former case, the coupling respects the internal dynamical structure of the 2LS and allows it to keep its 
coherence. 

To see how the distribution of modes in a field changes the picture, the cavity field calculation in Sec. 4 is useful. 
As shown in Fig. 1, the relaxation constant develops peaks and minima. The resonance effect is enhanced by a 
cavity size commensurate with the natural frequency of the 2LA and dissipation weakens. Narrow band resonance 
fluorescence as well as inhibition of spontaneous decay by frequent measurements - the Quantum Zeno effect - are 
interesting phenomena which our equations can provide finer details. 

Non-Markovian processes involve memory effects (nonlocal in time) . For the QBM problem, except for the case of 
high temperature Ohmic bath which gives Markovian dynamics, other types of spectral density (supraohmic) or at low 
temperatures, the dynamics of the system is non-Markovian [Q. When the reaction time of the bath is comparable to 
or faster than the natural time scale of the system (luq), one also expects to see non-Markovian behavior. By contrast, 
the 2LA is quite different: At zero temperature there is only one timescale = A _2 cj _1 >> lj^ 1 that determines 
both decoherence and relaxation. There is no memory effect and hence the process is Markovian. We expect that in 
finite temperature the dynamics of the 2LA will be nonMarkovian |3S[| , This is because there are more ways for the 
atom and the field to get entangled, and the memory effects of their interaction would presumably persist. 

In conclusion we find that the 2LS interacting with an EM field is far more coherent than what is commonly be- 
lieved, the misconception probably arising from the mistaken identification of this system with the Brownian model 
of an oscillator interacting with a harmonic oscillator bath. 
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APPENDIX A: ATOM-FIELD INTERACTION: TWO-LEVEL SYSTEM 



In this Appendix we give a rather detailed derivation of the Hamiltonian for a nonrelativistic atom interacting with 
a second-quantized electromagnetic field under the dipole, rotating wave and two level approximations. This is to 
facilitate the comparison of our model (II. 1) with a± coupling with that used by others (II. 3) with <r z coupling. (See 
Introduction). To make this also useful for later papers in this series, we have included atomic motion. Note that 
the convention here is closer to jjj than that used in the text which is closer to 0. The conversion is explained in 
footnote 2. 

The dynamics of a moving atom (mass M, momentum P) whose electrons (charge e, mass m) interact with an 
electromagnetic field (vector potential A, Coulomb potential V) is decribed by the (classical) Hamiltonian 

H=^ T + ^-(p-eA) 2 + eV(x)+H b (A.l) 
2M 2m 

where Ht is the Hamiltonian for the electromagnetic field. Expanding out one can write this as 

H = H a + H e + H b + H c (A.2) 

where 



describes nonrelativistic atom motion, 



2m 

decribes the dynamics of the (one) electron, while 



p2 

H a = (A.3) 



H e = ^- + eV( X ) (A.4) 



H D i = -— A p (A.5) 
m 



and 



H c2 = A 2 (A.6) 

describe the coupling between the electron and the field. The second term makes no contribution to one-photon 
processes and will be ignored. We will refer to Hq = H e + Hj, as the unperturbed Hamiltonian, and Hj = H c \ the 
interaction Hamiltonian. 

In a second-quantized form, the Hamiltonian for the radiation field is given by 

H b = ^kS k S k (A.7) 

k 

where 6k are the creation and annihilation operators for the kth normal mode of a free massless vector field. Thus 

for the field vacuum 6k|0) = 0, &k|0) = 0, [&k, tfj] = &c,k'i f° r au " k- We can perform a harmonic decomposition of the 
vector potential of the electromagnetic field 

i 

A(x, t) = J2 2 [ b ^ a (x)e^ 4 + 6 kCT Uk CT (x)e-^ 4 ] (A.8) 

where, assuming the field is contained in a box of size L, the spatial mode functions Uk CT is given by 

u kCT (x) = ir 3 / 2 e kCT ./k(x) (A.9) 
Here ek CT is the unit polarization vector and a = 1,2 are the two (transverse) polarizations. In free space, /k(x) = 

g-ikx 

We assume that electron motion is much faster than the motion of the atom, thus it sees a stationary central 
Coulomb potential around the center of mass of the atom. Denoting the (time- independent, nonrelativistic) electronic 
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wave function eigenstate belonging to the eigenvalue Ei by (pi, i.e, H e <pi = Ei<pi, we can write the Hamiltonian for the 
electron in the second-quantized form as 

H e = Y J E l a\a l (A. 10) 

i 

where i labels the bound states of the eletron (we assume vanishing probability for the atom to ionize) and a] and &i 
are the creation and annihilation operators. As fermions they obey the anticommutation relations d\dj + djd\ — 

To perform perturbation theory, the electronic wave function of the interacting system is expanded in terms of the 
eigenfunctions (pi of the unperturbed Hamiltonian, with basis formed by the direct product of the electron and the 
field states Thus the electron field operator V>(x) can be expanded as 

^)=^4i^(x). (A.ll) 

i 

With this, the interaction Hamiltonians is given by 

Hi = ~J ^ f (x)(p- A)i>(x)d 3 x (A.12) 



or in terms of dj , 6k operators 



Hi = h^2 a l % (flyk&k + 9ijkb£) (A. 13) 



where 



flyk 



1= [ <^(x)u kCT (x) • p^(x)d 3 x. (A.14) 
t^keo J 



Dipole Approximation 

Now consider conditions when the spatial variation of the vector potential A of the electromagnetic field is small 
compared to the electronic wave function ip, one can expand /k(x) in Uk CT (x) around the position of the atom 
x = X + <5x: 

e *k.x = gik-xjj + k . 5x _ I( k . 5x )2 + (A. 15) 

The dipole approximation amounts to keeping just the leading term. Doing so, we can take the field mode function 
/k(x) out of the integration above and evaluate it at the atomic position. To evaluate 



piP<j>jd 3 x 

m ' 



we make use of 



^ = ^i = l[^£ e ] (A.16) 
at m in 

yielding, 

— / (f>ip(j)jd 3 x = iuj^pdij (A. 17) 

TYl J 

where Uuiij — Ei — Ej, and dij = ej (piX(pjd 3 x is the dipole matrix element, = dji. Define 

d -^~V&V d --^ (A ' 18) 

Note that J^k = rfjik- With this under the dipole approximation, 
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9ij\n = rfyk/k(X) 



(A.19) 



Rotating Wave Approximation 

In the interaction picture, recalling that the time evolution of the ladder operators are given by 

a\{t) = a\e iuJit a^t) = aje'^ and b k (t) = b^e" 1 ^ 1 (A.20) 
the interaction Hamiltonian Hj in the interaction picture becomes 

&! = h 9ijka\ a b k e^-^ + h ]T g ij1t a\ a, b[ e ^+^ (A.21) 

i,j,k ijV. 

where We see there are two types of oscillatory terms present: e l ( w «± w k). Processes most effective 

in the absorption or emission of a photon by the atom correspond to those with near resonance frequency Ujj ~ c^k- 
Assuming u>tj > (Ei > Ej), the first type with e - J ( w ^+ w k) j^g a rapidly oscillating phase and its contribution is 
small compared with the second type with e - l ("ij-^k) wnose stationary phase at near resonance gives a large contri- 
bution. Physically the first type corresponds to either the excitation of the atom along with the emission of a photon 
or the relaxation of the atom along with the absorption of a photon, which is less probable than the second type 
corresponding to the excitation of an atom upon the absorpsion of a photon or the relaxation of an atom with the 
emission of a photon. We shall therefore ignore the first type of terms, which amounts to working under the Rotating 
Wave Approximation. This is the second major approximation in this standard model. 



Two-Level Atom 



Let us now consider the idealized case when the atom has only two electronic states, |+ >, |— > corresponding 
to i = 2,1, with energies equal to E± — ±i?kj - (The two states can interchangeably be labeled as |1 >, |0 >, or 
|e >, \g > or | |>, | !>•) Thus w, f2 i = u>o- Thus 

H a = -^-(a^ - a[ai) = —^~<Tz = nuj S z (A.22) 

where we have introduced a Pauli matrix (2x2) representation a z — diag(l, — 1) For the interaction Hamiltonian 
above, under the RWA, in the i = 2, j = 1 contribution to the summation, the first line containing e l ( UJ o-^ k ) \ s kept^ 
while the second line is dropped. The reverse is true for the i = l,j = 2 term. The interaction Hamiltonian (in the 
interaction picture) now becomes 

Hi = hJ2 [.92ik4ai& k e-^ k ^ 2l)t + $12^^-^+^] (A.23) 

k 

Introducing the Pauli matrix representation for the fcrmion operators a\a\ — > a + = S+ and a\a2 — > cr_ = S-, and 
defining g^ = d 2 ik = ^i2k (recall g^ = c?y7c/k(X)) we can write the interaction Hamiltonian (in the Heisenberg 
picture) in a simple form: 

= (s+& k / k (x) + s_&t/ k (x))) (A 24) 

k 

Therefore the total Hamiltonian for our model of a moving atom interacting with a quantum electromagnetic field 
under the dipole, rotating wave and two-level approximation is given by 

6 = S + hL °^ z + + 3k (A( X )^+ 6k + /k( x )^- & k)] ( A - 25 ) 
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APPENDIX B: AN OPERATOR PROOF OF THE MASTER EQUATION 



We use the resolvent decomposition of the propagator 

-iEt 



e -iHt 



f dEe~ %tjt 



Then by writing H = H a + Hj we can expand the resolvent and get 

(E - H)- 1 = (E — Ho)- 1 ^ — (E — Ho)- 1 !!!)- 1 
= (E — H )-\l + (E- Ho)- 1 !!! + ^{E - H^H^E - H^Hj + . . .) (B.2) 

When we act the expanded resolvent in any vector |0, i) = |0) ® \i), denotes an eigenstate of the Hamiltonian of 
the two level system) we see that 

only expanded terms that contain alternating sequences of cr + and er_ survive. This makes the summation much 
easier. Also if we note that 



(E - Ho)- 1 Y J 9kb k (E - Hor'gk^bU ) 

t 

9k | 



{E uio)- 1 TF^T— 1°> : = ( E ^r'FiE) (B.3) 



fe 

Hence we can compute the matrix elements by resumming the expansion 

(z,0\(E - H)- 1 ^) = E- 1 (B.4) 

(z,l\(E-H)- 1 \0,0)=0 (B.5) 

(z, 1\(E - H)- 1 ^, l) = (E-u - F(E))- 1 (B.6) 

MB-V-'M- ^^ X^* (b.7) 

The reduced density matrix propagator in the energy basis of the two- level atom and for the vacuum initial state 

is 

J(i,j;t\m,n;0) = J dEdE'e- l{E - E,)t J DzDze-^ k SkZk (B.8) 

{(z,i\(E - i/r^rnXO^lOE - Hy^zJ)) 

We can then verify that the only nonzero elements mn — > ij are the following and their conjugates 

J(0, 0; t|0, 0; 0) -» E^E'- 1 (B.9) 

J(0, 1; t|0, 1; 0) -» E~ 1 {E' - u - F(E'))- 1 (B.10) 

m o; ,10,0,0) ~ E _ E -° F(E) g * :? (g) Y. (E .J e . u) < B -"» 

J(ll;i|l, 1; 0) - (E - ^ - F(E))~ 1 (E -u - F(E'))- 1 (B.12) 

Then it is easy to check that this reproduces the propagation as given by equation (3.21) that was obtained through 
the influence functional method. Indeed 

is exactly the same as the one defined by equation (3.9). 
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